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The modified HFB (MHFB) theory at finite temperature is derived, which conserves the unitarity 
relation of the particle-density matrix. This is achieved by constructing a modified quasiparticle- 
density matrix, where the fluctuation of the quasiparticle number is microscopically built in. This 
matrix can be directly obtained from the usual quasiparticle-density matrix by applying the sec- 
ondary Bogoliubov transformation, which includes the quasiparticle occupation number. It is shown 
that, in the limit of constant pairing parameter, the MHFB theory yields the previously obtained 
modified BCS (MBCS) equations. It is also proved that the modified quasiparticle RPA, which is 
based on the MBCS quasiparticle excitations, conserves the Ikeda sum rule. The numerical calcula- 
tions of the pairing gap, heat capacity, level density, and level density parameter within the MBCS 
theory are carried out for 120 Sn. The results show that the superfluid - normal phase transition is 
completely washed out. The applicability of the MBCS up to a temperature as high as T ~ 5 MeV 
is analyzed in detail. 

PACS numbers: PACS numbers: 21.10.Ma, 21.60.Jz, 24.60.-k, 27.60.+J 



I. INTRODUCTION 



The finite-temperature Hartree-Fock-Bogoliubov (FT-HFB) theory has been successfully applied 
to highly excited nuclei Q, . It offers a fully self-consistent treatment of the interplay between 
single-particle, pairing as well as rotational degrees of freedom for nuclei in thermal equilibrium. 

A major drawback of this theory is the omission of fluctuation effects, which can be classified as 
quantal and statistical fluctuations. Quantal fluctuations arise from the mean-field approximation 
to the exact density operator V. As a result, the Hartree-Fock-Bogoliubov (HFB) density operator 
2?hfb violates the symmetries of the single-particle Hamiltonian H such as the conservation of parti- 
cle number and spin. However, quantal fluctuations decrease as the temperature increases. Various 
methods such as the Lipkin-Nogami 0] method, particle- number projection Q, angular-momentum 
projection |fj, particle-number conserving pairing correlations |fj, etc. have been proposed to 
eliminate quantal fluctuations. 

On the contrary, statistical fluctuations, which appear at finite temperature (T ^ 0), increase 
with increasing T |7], |g, |9j . Even the knowledge of the exact density operator T> would not eliminate 
statistical fluctuations from the FT-HFB theory. The omission of statistical fluctuation effect leads to 
the violation of another symmetry, namely the unitarity relation of the particle-density matrix Q,EJ. 
An immediate consequence of this symmetry violation is the collapse of the pairing gap A(T) at 
a critical temperature T c w ^A(T = 0) in all calculations for realistic nuclei within the FT-HFB 
theory and its limit, the finite-temperature BCS (FT-BCS) theory [H 13. Wi lie) . Such collapse of the 
pairing gap has been usually speculated as the signature of the superfuid-normal phase transition 
in finite nuclei. However, by using the Landau macroscopic theory of phase transition [ijj, Moretto 
has shown a long time ago that statistical fluctuations wash out such phase transition in finite 
systems such as nuclei, where these fluctuations are indeed quite large [12 • This conclusion has 
been confirmed recently by the calculations within the modified BCS (MBCS) theory 0,0]. The 
latter employs the modified quasiparticles obtained by a secondary Bogoliubov transformation of 
usual quasiparticles explicitly involving the quasiparticle occupation numbers. Other approaches 
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such as the static-path approximation |15|, llq . shell-model Monte-Carlo approach [13, modern 
nuclear shell model calculations 0, as well as the exact solution of the pairing problem also 
show that pairing correlations do not abruptly disappear at T 7^ 0. 

Another example of symmetry violation caused by the HFB and/or BCS theories is the violation of 
the Ikeda sum rule within the renormalized quasiparticle random-phase approximation (renormalized 
QRPA). The Ikeda sum rule states that the difference S~ — S + = (2J + l)(N— Z) between the total 
strength S~ of (3~ transitions and of (3 + ones, S + , is independent of models, where N and Z are the 
neutron and proton numbers, respectively, and J is the angular momentum of the transitions |20|. 
The renormalized RPA (or the renormalized QRPA, which includes pairing correlations), is an 
approach to take into account the Pauli principle between the particle (quasiparticle) pairs, which 
the RPA (QRPA) ignores This renormalises the RPA forward-going X and backward- 

going Y amplitudes as well as the two-body interaction matrix elements by a factor, which involves 
the particle (quasiparticle) occupation numbers in the correlated ground state. As a result, the 
collapse of the RPA (QRPA) at a critical value of the interaction parameter is avoided. However, 
it has been soon realized that the renormalized QRPA violates the Ikeda sum rule |24|. Several 
approaches were proposed recently to resolve this problem [2^, |2(| . 

The goal of this paper is to derive a modified HFB (MHFB) theory at finite temperature, which 
conserves the unitarity relation of the particle-density matrix. It will be shown that this can be 
achieved by using a modified quasiparticle-density matrix, which takes into account the statistical 
fluctuation of quasiparticle number microscopically. This modified quasiparticle-density matrix can 
be alternatively obtained by applying the secondary Bogoliubov transformation in Refs. [T3L H^j on 
the particle-density matrix at zero temperature. It will be demonstrated that the BCS limit of the 
MHFB equations yields the modified BCS (MBCS) equations, which have been obtained previously 
in Refs. [HQ- It will also be proved that the modified QRPA [3, obtained by using the MBCS 
quasiparticles, conserves the Ikeda sum rule. 

The paper is organized as follows. Sec. II summarizes the main features of the FT-HFB theory 
and its violation of the unitarity relation. The MHFB theory, which restores the unitarity relation, 
is derived in Sec. III. The MBCS equations are derived as the limit of the MHFB ones in the 
same section. The restoration of the Ikeda sum rule within the modified QRPA is shown in the 
Appendix. The theory is illustrated in Sec. IV by numerical calculations of the pairing gap and 
thermodynamic quantities such as the heat capacity, level-density parameter, and level density as 
functions of temperature for 120 Sn. The same section also discusses in detail the applicability of the 
MBCS equations in numerical calculations using realistic single-particle energies at high temperature. 
The paper is summarized in the last section, where conclusions are drawn. 



This section summarizes the main features of the FT-HFB theory, which has been derived by 
Goodman in Ref. [l| . They are essential for deriving the MHFB theory at finite temperature in the 
present paper. 



The HFB theory is based on the self-consistent Hartree-Fock (HF) Hamiltonian with two-body 
interaction 



where .. denote the quantum numbers characterizing the single-particle orbitals, %j are the 
kinetic energies, and Vijki are antisymmetrized matrix elements of the two-body interaction. The 
HFB theory approximates Hamiltonian by an independent-quasiparticle Hamiltonian -£/hfb 



II. REVIEW OF FT-HFB THEORY 



A. 



HFB Hamiltonian 




(1) 




(2) 
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where N is the particle- number operator, /i is the chemical potential, Eq is the energy of the ground- 
state |0), which is defined as the vacuum of quasiparticles: 

Oi|0)=0, (3) 

and Ei are quasiparticle energies. The quasiparticle creation a\ and destruction on operators are 
obtained from the single-particle operators a\ and <Zj by the Bogoliubov transformation, whose 
matrix form is 



U V 

v* u* 



(4) 



with the properties 



UU* + VV^ = 1 , UV T + VU T = , (5) 



where 1 is the unit matrix, and the superscript T denotes the transposing operation. The two- 
body interaction term of Hamiltonian Q, expressed in terms of quasiparticle operators using the 
transfromation Q, contains also the terms ~ a\a^a\a\ , ajajaj.az, a|aja;afc, and their hermitian 
conjugated parts. These terms are neglected in the HFB approximation. They play the role of 
residual interaction beyond the quasiparticle mean field. The quasiparticle energies Ei and matrices 
U and V are determined as the solutions of the HFB equations, which are usually derived by applying 
either the variational principle of Ritz or the Wick's theorem |5j. 



B. Thermodynamic and statistical quantities within FT-HFB theory 

At finite temperature T the condition for a system to be in thermal equilibrium requires the 
minimum of its grand potential il 

Q = S -TS - fj,N , (6) 
with the total energy £ , the entropy S, and particle number N, namely 

sn = o . (7) 

This variation defines the density operator T> with the trace equal to 1 

TrD = 1 , 5n/dV = (8) 

in the form 

V = Z -i e -P(H-»N) ) Z = Tr[e-^ H -^}, P = T-\ (9) 

where Z is the grand partition function. The expectation value -< O >- of any operator O is then 
given as the average in the grand canonical ensemble 

-< 6 y= Tr{v6) . (io) 

This defines the total energy £, entropy S, and particle number N as 

S=~<Hy= Tr(VH) , S= - ~< V\nV y= -Tr(V\nV) , N =~< N y = Tr{VN) . (11) 

The FT-HFB theory replaces the unknown exact density operator V in Eq. © with the approxi- 
mated one, 2?hfb, which is found in Ref. [l| by substituting Eq. J3J) in to Eq. © as 

2?hfb = HiriiAfi + (1 - -ft)] , (12) 
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where A/i is the operator of quasiparticlc number on the i-th orbital 

Mi = a\ ai , (13) 

and rii is the quasiparticle occupation number. Within the FT-HFB theory m is defined according 
to Eq. (JTUJl as 

* = <#<> = ^wti ' (14) 

where the symbol (...) denotes the average similar to IjlOf) . but in which the approximated density 
operator 2?hfb replaces the exact one, i.e. 

(6) = Tr(2?HFB0) • (15) 

That the quasiparticle occupation number m at finite temperature is given by the Fermi-Dirac 
ditsribution as in Eq. I|14l) within the framework of the independent-quasiparticle approximation 
has been also proved a long time ago by Zubarev using the double-time Green function method [23 
(See also the Appendix A of Ref. [14j)- The quasiparticle energy Ei in Eq. I|14l) is found by solving 
the FT-HFB equations summarized in the next section. 



C. FT-HFB equations 



The generalized particle-density matrix R is related to the generalized quasiparticle-density matrix 
Q through the Bogoliubov transformation Q as 

R=U^QU , (16) 

where 

«=(4- !-',)-(; i-„). m 

with 

/ rr* tv - * \ 

= ( y V v ) , Utf = l. (18) 

The matrix elements of the single-particle matrix p and particle pairing tensor r within the FT-HFB 
approximation are evaluated as 

Pij = (a]ai) , = (aj-ffli) , (19) 

while those of the quasiparticle matrix q are given in terms of the quasiparticle occupation number 
since 

qij = (a'jOi) = SijUi , = (ajOti) = , (20) 

which follow from the HFB approximation J2J. Using the inverse transformation of the particle 
densities are obtained as 

p=U T nU* + V^(l-n)V , r = U T nV* + V f (l - n)U . (21) 

By minimizing the grand potential Q according to Eq. JJJ, Goodman has derived in Ref. the 
FT-HFB equations in the following form 

H A \(U i \ =E (UA > (22) 



A* -W* M Vi / l \V t 
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where 

H = T - fi + T , Tjj = Vikjipik , A,j = i VijkiTki ■ (23) 



2 

fcZ fci 



(For the details of the derivation see Sec. 4 of Ref. The total energy £, entropy S, and particle 
number N from Eq. are now given within the FT-HFB theory as 

£ =Tr[(T+ir)p+^Art] , (24) 



5 = - ^Klnn, + (1 - ni)ln(l - m)} , (25) 



N = Trp , (26) 

from which one can easily calculate the grand potential fi (JBJ . 
In the limit 



= -G« , (27) 



i ij j l 3 ' 



where |«) denotes the time-reversal state of \i), Eqs. and iffiH) yield the well-known FT- 

BCS equations. For spherical nuclei and with all the pairing matrix elements equal to Gy = G, the 
FT-BCS equations have the form 

A = G ^i u 3 v 3 (! - 2n j) > ( 28 ) 



AT = 2^%[(l-2n>? + n J -] , (29) 
i 

where 2f2j = 2 j + 1 is the shell degeneracy. The quasiparticle energies , and Uj and Uj coefficients 
are given as 



D. Violation of unitarity relation within FT-HFB theory 

At zero temperature (T = 0) the quasiparticle occupation number vanishes: =0, and the average 
l|15(l reduces to the average in the quasiparticle vacuum The quasiparticle-density matrix Q Ill7|l 
becomes 

Q(T = 0) = Q = ( ° J ^ , for which Q 2 = Q ■ (31) 

Therefore, for the generalized particle-density matrix Rq = R(T — 0) the following unitarity relation 
holds 

Rq = Ro , (32) 

where 

R Q =U^Q U. (33) 
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However, the idempotent 132f) no longer holds at T ^ 0. Indeed, from Eqs. (|16fl and l|17|) it follows 
that 

R- R 2 =U\Q -Q 2 )U , (34) 

which leads to 

Tr(i?-i? 2 ) = Tr(Q-Q 2 ) = 2^n 4 (l-n 2 ) = 2(<5A0 2 ^0, (T ^ 0) . (35) 

i 

The quantity <W 2 = ~ n i) m Eq. I|35|l is nothing but the quasiparticle- number fluctuation. 

This can be easily checked by calculating 

sm 2 = (m 2 ) - (ao 2 = (Em + J2ttA) - XX - £ n ^ 

= XIn i (l-n i )=X;^, (36) 

i i 

where 

5M 2 =n l {l-n l ) (37) 

is the fluctuation of quasiparticle number on the i-th orbital. We've just seen that the violation 
of the unitarity relation (|32|) for the generalized single-particle density matrix R occurs at T ^ 
due to the fact that the HFB approximation @ and the density operator 2?hfb H2|) exclude the 
quasiparticle- number fluctuation (|36|l from the quasiparticle-density matrix l|17|l [9(- Therefore, in 
order to restore the idempotent of type (|32|l at T =^ a new approximation should be found such 
that it includes the quasiparticle-number fluctuation [Eqs. (|36|) and l)37|) ] in the quasiparticle-density 
matrix. 

III. MODIFIED HFB (MHFB) THEORY AT FINITE TEMPERATURE 

By including the quasiparticle-number fluctuation (|36l) . a part of the higher-order terms ~ 
alajCtkCti, neglected as the residual interaction beyond the FT-HFB quasiparticle mean field, will 
be taken into account. As the result the mean field of usual quasiparticles itself will be modified. 
This leads to the new quasiparticle energy Ei and chemical potential ft, which will be found as the 
solution of the modified HFB (MHFB) equations to be derived in this section. 

A. Restoration of unitarity relation 

Let us consider, instead of the FT-HFB density operator 2?hfb l|12() . an improved approximation, 
T>, to the density operator T>. This approximated density operator T> should satisfy two following 
requirements: 

(al) The average 

((d)) = Tr(v6), (38) 
in which V is used in place of T> (or 2?hfb ) , yields 

R=lfiQU (39) 
for the Bogoliubov transformation U (|18() , where one has the modified matrices 

*=(4*i-V)' M-pi-V)' (40) 
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with 



Pij 



((OjOi)) 



(41) 



Qij = ((<4a*)) = Si- 



Uj = {(ajOi)) 



(42) 



instead of matrices R and Q in Eqs. (|17|l . I|19|) . and (|2U|) . The non-zero values of t%j in Eq. i|42|l are 

caused by the quasiparticle correlations in the thermal equilibrium, which are now included in the 
average ((■•■)) using the density operator V. 

(a2) The modified quasiparticle-density matrix Q satisfies the unitarity relation 



(<2) 2 = Q ■ 

The solution of Eq. (|43|l immediately yields the matrix A in the canonical form 



(43) 



-Ai\ 
Ai 





A 2 



-A 2 




A; 



\/ Ui(l - Tli) 



(44) 



V / 

Comparing this result with Eqs. 1)36(1 and (|37[1 . it is clear that tensor A consists of the quasiparticle- 
number fluctuation <Wi = yjfii(l — fii). From Eq. I|39|) it is easy to see that the unitarity relation 
holds for the modified generalized single-particle density matrix R since 



R-R 2 =W f (Q -Q 2 )U = 



(45) 



due to Eq. I|43|) and the unitary matrix hi. 

Let us define the modified-quasiparticle operators a\ and on, which behave in the average l|38(l 
exactly as the usual quasiparticle operators ot\ and on do in the quasiparticle ground state, namely 



{{a\a k )) = «q!4)) = ((Sfcfii)) = . 



(46) 



In the same way as for the usual Bogoliubov transformation J3J, we search for a transformation 
between these modified-quasiparticle operators (aj, on) and the usual quasiparticle ones (a\, ai) in 
the following form 



o 



W 2 

z w 



a 
a 



with the unitary property similar to Eq. |[5jl for U and V" matrices 

wai' + zz^ = 1 . 

Using the inverse transformation of (|47|l and the requirement we obtain 

rii = {{a\ai)) = ^ 



c ik ■ 



(47) 



(48) 



(49) 



From this equation and the unitary condition l|48|) . it follows that zz* = n and ww^ = 1 — n. Since 
1 — h and n are real diagonal matrices, the canonical form of matrices w and z is found as 



(wi 
wi 



VI -n = 



w 2 
w 2 





-21 



-* 2 




(50a) 
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where 

Wi = y/l - Hi , Z l = yWi . (50b) 

We see now that, just like Eq. Q}, which is the generalized form of the Bogoliubov transformation 
for the BCS case, Eq. (|47|l with matrices w and z uniquely defined in Eq. I|5(J|I is the genearlized form 
of the secondary Bogoliubov transformation used in Refs. [l3l[T3|. It expresses a simple relationship 
between the modified-quasiparticle operators (aj, cti) and the usual-quasiparticle operators (a\, a<) 
in the same fashion as that between the latter and the single-particle operators in the Bogoliubov 
transformation (0J. 

We now show that we can obtain the idempotent R 2 = R 1(45(1 by applying the secondary Bo- 
goliubov transformation 1)47(1 . which automatically leads to Eq. ((43|) . Indeed, using the inverse 
transformation of l(47() with matrices w and z given in Eq. (|50fl . we found that the modified 
quasiparticle-density matrix Q can be obtained as 



where 



and 



W^Q W=( * "_" )]t ) S g, (51) 

\ \/n{l — n) 1 — n J 



«ata)) ((ad)) \ _ / \ = 2 _ = 
((at a t)) 1 - «ata)) J ~ ^ 1 J ' y ° " ^° ' 

due to Eq. 1(46(1. This result shows another way of deriving the modified quasiparticle-density 
matrix Q 1)40(1 from the density matrix Qq of the modified quasiparticles (aj, Si). This matrix Qo is 
identical to the zero-temperature quasiparticle-density matrix Qq 1(31(1 . Substituting this result into 
the right-hand side (r.h.s) of Eq. 1(39(1 . we obtain 

R = WQ U , (54) 

where 

This equation is the generalized form of the modified Bogoliubov coefficients Uj and Vj given in 
Eq. (6) of Ref. [3 or Eq. (38) of Ref. 0. From Eqs. (JTHJ, and JS^Jl, it follows that 

UU^ = 1, i.e. transformation ((54(1 is unitary. Therefore, from the idempotent ((53(1 it follows that 
R 2 = R. We have just shown that the secondary Bogoliubov transformation ((47(1 allows us to 
take into account the fluctuation of quasiparticle number and restore the unitarity relation of the 
generalized particle-density matrix . In this sense, the approximation discussed in the present 
section is a step beyond the thermal mean field of usual quasiparticles. As the result, the thermal 
quasiparticle mean field, which was defined within the FT-HFB approximation, is modified due to 
thermal quasiparticle-number fluctuation. 2 



1 An alternative approach to the unitarity problem was propposed in Ref. l2Sil making use of the thermo field 
dynamics l2fl . 

2 An exact theory on quasiparticle excitations at T = should define the vacuum and quasiparticles in terms of 
exact eigenstates of the many-body system But in this case, a simple mathematical relationship between the 
exact quasiparticles and the usual particles of the system no longer exists. The advantage of the Bogoliubov- 
type quasiparticles is the linear relationship between them and the usual particles. However, the corresponding 
vaccum and single-quasiparticle state are now only approximations of the exact eigen functions of the many-body 
Hamiltonian. Similarly, at T ^ 0, when the average over the individual compound systems is replaced by that over 
the grand canonical ensemble I1UI . the density operators I^hfb 1121 and T> 1381 are different approximations of the 
exact density operator T> 
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B. MHFB equations at finite temperature 

With all the thermal degrees of freedom now included in U, Eq. 1(54(1 formally looks the same 
as the usual HFB approximation at T = 1)33(1 . which connects Ro to Qq. Applying the Wick's 
theorem for the ensemble average |5|, one obtains the expressions for the modified total energy £ 



£ = Tr[(T+-f)p- 



;Aft] , 



where 



^ Vikjipik 

kl 



(56) 



(57) 



It. 



Vijkl T kl 



(58) 



kl 



From Eq. (|54|l we obtain the modified single-particle denity matrix p and modified particle-pairing 
tensor f in the following form 



p = U T nU* + V f (l - n)V + U 1 



\J n(l - n) 



t 



V + V^n(l - n)U* , 



(59) 



U T nV* +V\l-n)U + U 1 



1 1 



U + V^ y / h(l~n)V* 



(60) 



As compared to Eq 
two terms 



(|21|) within the FT-HFB approximation, Eqs. 
It 



l|59|l and l|6U|) contain the last 
[a/ 77,(1 — fi)y and ~ \J n(l — fi), which arise due to quasiparticle-number fluctuation. 
Also the quasiparticle occupation number is now n [See Eq. H42|l ] instead of n (|14|l . 

We derive the MHFB equations following the same variational procedure, which was used to derive 
the FT-HFB equations in Sec. 4 of Ref. According it, we minimize the grand potential SCl = 
by varying U, V , and n, where 



n = £ -TS-pN 



(61) 



Due to Eq. (JSJ), the variations SU and SV are not independent. They are found by using an 
infinitesimal unitary tranformation of The obtained infinitesimal variations U' = U + 8U and 
V' = V + SV together with n' — n + Sn are then used in Eqs. (|59|l and (|60|l to obtain p = p + Sp 
and f'=f+ St. Substituting them into Eq. Q61[l one obtains f2' = f2 + SCl, where <5f2 is expressed 
in terms of Sp, St, and Sn as independent variations. By requiring the coefficients of Sp and Sf 
vanish and following the rest of the derivation as for the zero-temperature case, we finally obtain 
the MHFB equations, which formally look like the FT-HFB ones 



H A 

-A* -H* 



Ui 



= Ei 



(62) 



where, however 



H = T-fl 



(63) 



with r and A given by Eqs. (|57(l and (|58(l . respectively. The equation for particle number N within 
the MHFB theory is 



N = Trp 



(64) 
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By solving Eq. I|62l) . one obtains the modified quasiparticle energy Ei, which is different from Ei 
in Eqs (1221) and/or (|30|l due to the change of the HF and pairing potentials. Hence, the MHFB 
quasiparticle Hamiltonian Hmhfb can be written as 

H - fiN » i/MHFB = E Q + Eitii , (65) 

i 

instead of © . This implies that the approximated density operator D (|38|) within the MHFB theory 
can be represented in the form similar to 1|12[1 . namely 

V = Pmhfb = Y[[niATi + (1 - -&)]■ (66) 

i 

From here it follows that the formal expression for the modified entropy S is the same as that given 
in Eq. J25J|, i.e. 

S = -Qfiilnfti + (1 ~ " »<)] . (67) 

i 

Using the thermodynamic definition of temperature in terms of entropy T — SS/SS and carrying 
out the variation over dfii, we find 

8£ i=. „ SS „, / 1 — fii \ , . 

— = E i = T- r - =Tln • 68 

orij orii \ Ti-i y 

Inverting Eq. (|68|l . we obtain 

This result shows that the functional dependence of quasiparticle occupation number fii on quasipar- 
ticle energy and temperature within the MHFB theory is also given by the Fcrmi-Dirac distribution 
of noninteracting quasiparticles but with the modified energies Ei defined by the MHFB equations 
(|62H . Therefore in the rest of the paper we will omit the bar over fii and use the same Eq. I|14|) 
with Ei replaced with Ei for the MHFB equations. 



C. Modified BCS (MBCS) theory at finite temperature 

1. MBCS equations 

The MBCS equations at finite temperature have been derived previously in Refs. 0, Q] using 
the secondary Bogoliubov transformation l|47|l for the BCS case. We will show below that these 
MBCS equations emerge as the limit of the MHFB equations derived in the preceding section. 

In the BCS limit H27|) with equal pairing matrix elements Gij = G, neglecting the contribution of 
G to the HF potential so that f = 0, the HF Hamiltonian becomes 

Hy = (e< - H)S tj . (70) 



3 Note that, there remains the residual interaction, even due to pairing alone, beyond the MHFB quasiparticle 
mean fi eld . A t T = this can be treated as ground state correlations within the renormalized and/or modified 
QRPA [HEl- 

As the result n% deviates from the Fermi-Dirac distribution, especially if different multipolarities 
of the two-body residual interaction are taken into account. However, for the monopole pairing interaction alone 
as considered in this paper, such deviation is negligible [See Appendix B of Ref. 1 141 1. At T ^ the quasiparticle- 
number fluctuation beyon d th e quasiparticle mean field leads to the entropy effect within the renormalized RPA, 
which was studied in Ref. I.'iil . 
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The pairing potential l|58|) takes now the simple form 

a = -G5>*fc- (71) 

fc>0 

The Bogoliubov transformation f° r spherical nuclei reduces to 

a jm = U j a ]m + v j(~y + a j-rn , 

(-y +m a^ m = Uj (-y +m aj . m - Vj at m , (72) 
while the secondary Bogoliubov transformation (|47|l becomes 

& ]m = V 1 - n J a ]m - Vnj(-y +m aj-m , 
(-y+^-m = ^T^-(-y +m a^ m + Jn s a\ m . (73) 



The U, V, 1—n, n, and \J n{l — n) matrices are now block diagonal in each two-dimensional subspace 
spanned by the quasiparticle state \j) and its time-reversal partner \j) = (— ) J+m |j — to) 



1 — n 



1 - m o \ / 



1 — rij J ' V rij 



Substituting these matrices into the r.h.s of Eqs. I|59|l and (|6Uf) . we find 



(75) 



p. ~. = (1 - 2rij)«j + «j - ZyJrijil-n^UjVj , (76) 



-(1 - 2nj)ujVj + \Jnj(l - nj)(uj - uf) . (77) 



Substituting now Eqs. (|77|l and (|76[) into the r.h.s of Eqs. (|71ll and (|64|l . respectively, we obtain the 
MBCS equations for spherical nuclei in the following form: 



A = G^fij-Kl - 2n i )u i »j - J nj (l - - i, 2 )] , (78) 



jV = 2^0^(1 - 2nj)v] + n 3 - 2J nj {l - n^UjVj] . (79) 



Equations (|78|l and l|79l) are exactly the same as the MBCS Eqs. (23) and (24) in RefjRl or E q s - 
(39) and (40) in Ref. |14|. We've just shown that the MBCS equations in Refs. [Tll^ emerge as 
the natural limit of the MHFB equations at finite-temperature. 
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For convenience in further discussions we rewrite the MBCS gap in Eq. (|78|l as a sum of quantal 
Aq and thermal-fluctuation <5A parts as 

A = A Q + 6 A , (80) 

where the quantal gap Aq is 

A Q - ^(A Q ), , (A Q ), = GfljUjVjil - 2nj) . (81) 

3 

It is called quantal since it is caused by quantal effects starting from T = 0, where it is equal 
to the BCS gap, and decreases as T increases because the Pauli blocking becomes weaker. The 
thermal-fluctuation gap 5A, referred to hereafter as thermal gap, is given as 

5 A = SAj , 6Aj = Gttjiv? - u])8N , (82) 

3 

and arises due to the thermal quasiparticle-number fluctuation 8Mj at T ^ 0. Therefore, comparing 
the FT-BCS equations J2HJ an d J5*jJ with the MBCS ones, O and we see that the latter 

explicitly include the effect of quasiparticle-number fluctuation ~ SAfj <|37|) in the last terms at 
their r.h.s, which are the thermal gap (|82(l in Eq. i|78|) and the thermal-fluctuation of particle 
number SN = J2j SN 3 = i Y. J i --/".i r j(*Xi ] in E( l- CB- These terms 

are ignored within the 

FT-BCS theory. Hence Eqs. (|78l) and (|79|l show for the first time how the effect of statistical 
fluctuations is included in the MBCS (MHFB) theory at finite temperature on a microscopic ground. 
So far this effect was treated only within the framework of the macroscopic Landau theory of phase 
transition fl^ . 



2. Thermodynamics quantities 



The total energy £ is found as 



£ — 2 Qj e j [(1 — 2 n j) v j + 11 j ~~ 2W7ij(l — n,j)ujVj] 



A2 
G 



(83) 



The heat capacity C is calculated as the derivative of energy £ ijBBil with respect to temperature T 

d£ 



c = 



dT 



The level-density parameter a is defined by the Fermi-gas formula 

E* _ £{T) - £(0) 



T 2 



T 2 



(84) 



(85) 



where E* = £(T) — £(0) is the excitation energy of the system. The quasiparticle entropy (|67ll is 
written for spherical nuclei as 



S = -2^2 ^3 [ n 3 lnn 3 + (1 - nj)ln(l 



2^3 



0E j 



2 0E 3 + ! 



ln(l + e 



-f3Ej > 



(86) 



Using the MBCS equations (|T5jl and (JJSJ, Eqs. iJHUJl, and (j5B*jl together with the expressions for 
Ei, ui, and Vi, which are the same as in Eq. 1)30(1 (with Ei replacing Ei and A replacing A), we 
found that the formal expression for the grand potential $ is also the same as that given within the 
FT-BCS theory [U 112 , namely 



4> 



(87) 
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The level density p(N, Z) is calculated as the inverse Laplace transform of the grand partition 
function e* . It is approximated as [32, 



p(N,Z) 



2tt\/2^D ' 



(88) 



where S = S n + Sz is the total entropy of the system, D is the determinant of the second derivatives 
of the grand partition function taken at the saddle point. It is given as 



D 



<9 2 <I>. 

da 2 N 



D; 



<9 2 $. 

do? z 



D 



N 



D, 



da' 2 daid(3 



Jaid/3 83 2 



i = N,Z 



(89) 



The formal expressions for the derivatives in the determinant are the same as given in Eqs. 
(B.15) - (B.17) of Ref. j^]. However, the derivatives of the gap A entering in these expressions are 
more complicate due to Eq. (|78|l . They are obtained here as 



<9A 

da 



rj -i A 2 



8 E 2 



2a j , 



(90) 



dA 
~df3 



Ej ttA&Kej - fiK^aj - fib,) + A 2 a 3 ] - f A 2 ^ + (e, - fi)[(ej - /2) ej + A 2 



/? V . %[£ 2 ^ + A 2 K- - &,)] -AV. %( £j - p)(/3 Cj + 



(91) 



where 



sech 2 (z) tanh(z) sech(z)tanh(z) 1 

^ = 2£ 2 ' J = ~jEf~ ' ° j = 2E] ' Z = 2^ j • ( j 

We've just derived the MHFB theory at finite temperature, which includes the quasiparticle-number 
fluctuation to preserve the unitarity of the modified generalized particle-density matrix. We've 
shown that the limit of this MHFB theory reproduces the MBCS equations obtained previously in 
Ref. [H|ll. For the sake of completeness we give in the Appendix the proof that, by using the 
secondary Bogoliubov transformation {73J, the modified QRPA indeed conserves the Ikeda sum rule. 



IV. ANALYSIS OF NUMERICAL RESULTS 



As an illustration for the modified HFB theory at finite temperature, we now discuss in detail 
the results of numerical calculations within its limit, the MBCS theory, of the pairing gap, heat 
capacity, level-density parameter, and level density for 120 Sn. The single-particle energies 6j used in 
the calculations are obtained within the Woods-Saxon potential at T =0. These discrete neutron and 
proton spectra include not only bound but also quasi-bound levels, which span an energy interval 
from around -40 MeV up to around 17 MeV. They include all the major shells up to N(Z) = 126 
as well as several levels in the next major shell N{Z) — 126 - 184 up to lfci7/2 orbital. They are 
assumed here to be independent of T. This assumption is supported by the results of temperature- 
dependent HF calculations, which show that for T < 5 MeV the variation of the single-particle 
energies with T is negligible |3J|. The value G„ = 0.13 MeV is adopted for the neutron pairing 
parameter so that the gap A„ for neutrons is about 1.4 MeV at T = 0. 
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A. Temperature dependence of pairing gap 

1. Open-shell case: Neutron pairing in 120 Sn 

Since the modified gap A is a function of T, the last term 5 A Ij82(l at the r.h.s of Eq. (|78|l raises a 
question about the validity of MBCS equation at high temperature. In fact, at first glance, it seems 
that, if the single-particle spectrum is such that 6 A is negative and its absolute value is greater than 
that of the first term at the r.h.s of Eq. H78fl a t a certain value of T, the gap A will turn negative 
and the MBCS approximation breaks down. In this section we will show that this does not happen 
in numerical calculations using the entire single-particle energy spectrum. 

Shown in Fig. 1 are the Bogoliubov coefficients Uj, Vj, quasiparticle occupation number rij, 
together with the combinations UjVj, uj — vj, 1 — 2rij, and SJ\fj as functions of single-particle energy 
tj for neutrons at several temperatures. These quantities determine the behavior of the gap 1781) 
as a function of T. They are rather symmetric functions from both sides of the chemical potential 
ft. The latter varies weakly around -6 MeV as T increases. The product UjVj decreases quickly 
with increasing T. At T = 5 MeV, it remains effective only in the region of ±5 MeV around 
ft. The difference uj — vj, which enters in the thermal part 6 A, remains rather insensitive to the 
variation of T. In general, the effect of pairing on the Bogoliubov coefficients Uj and vj and their 
combinations, UjVj, and uj — vj, is significant only in the region of at most ±10 MeV around ft. This 
situation is rather similar to that obtained within the BCS theory. However, for the quasiparticle 
occupation number rij and its combinations, 1 — 2rij, and SAfj, the situation is different. Here, with 
increasing T, these quantities, although having a peak near ft, spread over the whole single-particle 
spectrum as shown in Fig. 1 (d) - (f). For the quantal component (Aq).,- the maximum of UjVj 
comes always with the minimum of (1 — 2n,j) near ft. Beyond this region the product UjVj(l — 2rij) 
is small. However, for the thermal-fluctuation part of the gap, both regions far above and below 
ft are important. This means that, in difference with the BCS theory, where one can restrict the 
calculations with valence nucleons on some closed-shell core by renormalizing the pairing parameter 
G„, the calculations for open-shell nuclei within the MBCS theory are necessary to be carried out 
using the entire single-particle spectrum. 

This observation is demonstrated on Fig. 2, where the partial quantal (Aq>)j and thermal 5Aj gaps 
are shown as functions of single-particle energy ej at several temperatures. The quantal part (Aq)^ 
is always larger around ft, but its magnitude quickly decreases as T increases. On the contrary, the 
thermal part 5Aj is positive at tj < ft, and negative at 6j > ft. Its absolute value sharply increases 
with increasing T. In a realistic spectrum the number of single-particle levels below ft is usually 
larger than that of those above it. In the present example of 120 Sn, within the same energy interval 
of ~ 20 MeV from ft, the one below ft has 12, while the one above ft has only 8 single-particle levels. 
Therefore, the sum of partial thermal gaps SAj has more components in the region below ft, where 
the difference vj — uj is positive. As the result, by summing over all single-particle levels weighted 
over the shell degeneracy fij, the ensuing thermal gap <5A <|S2l) is always positive. 

Shown in Fig. 3 are the quantal Aq (dashed line) and thermal SA (dash-dotted line) gaps together 
with the total MBCS gap A (thick solid line) as functions of T. The BCS gap is also shown as the 
dotted line for comparison. It collapses at a critical temperature T c ss 0.79 MeV. This value almost 
coincides with the temperature of superfluid - normal phase transition estimated for infinite systems, 
which is about 0.567A(T = 0) On the contrary, within the MBCS theory, the quantal gap Aq 
never collapses, but decreases monotonously with increasing T. The thermal component 8 A increases 
first with T at T < 1 MeV, then starts to decrease with increasing T further, but still does not vanish 
even at T ~ 5 - 6 MeV. As the result, the total MBCS gap A has a temperature-dependence similar 
to that of the quantal gap Aq, except for a low-temperature region 0.5 MeV < T < T c , where it 
increases slightly with T because of the thermal gap 8A. As high T, the total gap A decreases 
monotonously with increasing T. This yields a long tail extending up to T ~ 5 - 6 MeV. 

In order to see how the change of configuration space affects the calculation of the MBCS gap, 
we also carried out several tests using cut-off spectra. Examples are shown in Fig. 4. The dashed 
line is the neutron gap obtained in the MBCS calculation after removing three lowest major shells 
(up to N = 28) from the single-particle energy spectrum. The calculations are then carried out 
by putting N = 42 particles on the N =28 core. The balance in the sum over the single-particle 
levels is lost with less levels below ft participating in the summation. The symmetry of the spectrum 
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with respect to p is destroyed. The gap collapses again, but at a much higher temperature T w 4 
MeV, although up to T ~ 2.5 MeV its temperature dependence is almost the same as that obtained 
using the entire spectrum. Removing from the other side of p two highest levels I&17/2 and liix/s 
makes the reduced spectrum rather symmetric again with respect to p. The balance in summation 
of SAj is restored. As the result, the temperature dependence of the gap is recovered as shown by 
the thin solid line. However, if one removes further one more level, namely the lj'15/2 one, i.e. the 
reduced space consists of only three major shells, 28 - 50, 50 - 82, and 82 - 126, the balance is 
destroyed again with more weight toward the positive values of SAj. The reduced spectrum now 
spreads from around -17 MeV up to around 1.6 MeV, which is strongly asymmetric with respect to 
p. In consequence, the high-temperature tail of the gap becomes much more enhanced as shown by 
the dash-dotted line in Fig. 4. Other tests using a 8-neutron, 20-neutron, and 50-neutron cores also 
show a similar feature. In these tests the parameter G v is renormalized so as to obtain the same 
value for A„(T = 0). With such renormalization of G v the BCS gap always remains the same. 

These results show the difference in practical calculations within the BCS and MBCS theories. In 
the BCS case, the calculation of the gap using a closed-shell core with a simple renormalizition of the 
pairing parameter G v yields the same result as that obtained using the entire single-particle energy 
spectrum. In the MBCS case, the most reliable way is to use the entire or as larger as possible single- 
particle spectrum. If using a limited spectrum is unavoidable, care should be taken to maintain the 
balance in the summation of partial thermal gap SAj . Otherwise, a resulting collapse or an enhanced 
tail of the gap in the high-T region would be simply an artifact caused by a limited space. As a matter 
of fact, a criterion for a good reduction is that the cut-off spectrum should be rather symmetric with 
respect to the region where the quantal pairing correlations are strongest, namely from both sides of 
the chemical potential, so that the effect of quasiparticle-number fluctuation is properly taken into 
account (See Figs. 1 (d) - (f) and 2 (b)). It is worth noticing that the limitation of the configuration 
space also yields a wrong behavior of the specific heat. This effect is known as the Schottky anomaly, 
according to which the specific heat reaches a maximum at a certain temperature and decreases as 
temperature increases further |35j . 

2. Closed-shell case: Thermally induced pairing correlations for protons in 120 Sn 

The MBCS gap equation l|78[l also implies that, in principle, thermal fluctuations can induce 
pairing correlations even for closed-shell (CS) nuclei. However, the situation here is different from 
that of the open-shell nuclei because of a large shell gap between the highest occupied (hole) orbital 
and the lowest empty (particle) one, which is about 6 MeV for protons in 120 Sn. At T — all the 
orbitals below p are fully occupied (Vj h =l, Uj h — 0, 6j h — fx < 0), while those above p are empty 
{vj h =0, Uj h = 1, €j p — > 0). Therefore the quantal gap Aq i|%T|) is always zero. Pairing is so 
weak that no scattering into the next major shell (particle orbitals) is possible. In such situation 
the approximation of the same pairing matrix elements may not be extended across a too large shell 
gap separating hole and particle orbitals, especially when fj h = 1 — rij h ^> rij p = fj , where fj is 
the single-particle occupation number. This restricts the summation at the r.h.s of Eq. (|78|l to be 
carried out at most over only the hole states. The MBCS gap A in this case is solely determined by 
the thermal gap SA (I82fl due to the quasiparticle-number fluctuation, namely 

Acs - SA CS a G^n^yjnj^l - n Jh ). (93) 

The thermally induced gap A,r for closed-shell proton system (Z = 50) in 120 Sn, obtained using 
Eq. (|93|l with the same value of pairing parameter as that for neutrons, G^ = G„, is plotted as a 
function of temperature in Fig. 5. This figure clearly shows that the pairing gap for a closed-shell 
system is different from zero at T 7^ and increases as T increases. However, its magnitude, which 
reaches a value of only around 2.6xl0 -5 MeV at T = 5 MeV, is practically negligible as compared 
to A v . Therefore we will put A w equal to zero in further discussions. 
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3. Comparison between microscopic and macroscopic descriptions of thermal fluctuation 

The effect of thermal fluctuations on the pairing gap was first studied using the Landau macro- 
scopic theory of phase transition ^lj by Moretto in Ref. Within the Landau theory, $ 1)87(1 
is treated as a function of the independent parameter A. The probability that the nucleus has any 
given value of A for the pairing gap is determined by the isothermal distribution 

P(A) oc e*< A ) . (94) 

The averaged gap (A) is calculated as ^2 

f n °° AP(A)dA 

This approach does not include quantal fluctuations. Therefore, as has been pointed out in Refs. 
0, IllL Il2| , at very low temperature or if nonequilibrium states vary too rapidly with time, quantum 
fluctuations dominate and Eq. i|95|) is no longer meaningful. 

The probability distribution P(A) (|94f) . calculated using the same neutron single-particle spectra 
for 120 Sn and the same pairing parameter G v , is plotted as a function of A at low and high temper- 
atures in Fig. 6 (a) and (b), respectively. At very low temperature, the most probable value, which 
is the BCS gap, coincides with the averaged one, resulting in a Gaussian-like shape with a peak at 
the BCS value of A(T = 0) ~ 1.4 MeV. As T increases, the distribution becomes skewed toward the 
lower values of A. Its maximum, which still corresponds to the solution of the BCS equation, moves 
to lower A and reaches A = at T = T c . This is shown in Fig. 6 (a), which is very similar to what 
obtained before in Fig. 1 of Ref. [12| for a uniform spectrum. As T increases further, the maxi- 
mum of the distribution still remains at A = 0, while its width continues to increase, showing the 
increase of thermal fluctuations. At hypothetically high temperatures (Fig. 6 (b)) the distribution 
approaches a Gaussian one in the following form 



I / A 2 \ IGT 

P(A)| T _~^ex P (-— j with a=^—. (96) 



This behavior means that Eq. I|94|l assumes the effect of thermal fluctuations on the pairing gap to 
be chaotic at both low and high T, By substituting Eq. I|96|) into the r.h.s of Eq. Q95fl. the integrals 
can be carried out analytically. The result is 



<A)| T ^oo = ^. (97) 

This result reveals the increase of the averaged gap (A) with increasing T at very high T within the 
Landau theory using the probability function (|94|) . 

The temperature dependence of (A) for neutrons in 120 Sn is displayed in Fig. 7 in comparison 
with the MBCS gap A. It is seen from this figure that the agreement between the microscopic 
treatment of thermal fluctuations in the pairing gap within the MBCS theory and the macroscopic 
one can be called at best qualitative. The gap does not collapse in both treatments, but while the 
tail of MBCS gap A clearly decreases at high T with increasing T, the temperature dependence of 
the averaged gap (A) remains rather flat, and even starts to increase slightly with T already at T > 
1 MeV because of Eq. This yields a (A) of about 0.6 MeV even at T ~ 4 - 5 MeV. In the 

low temperature region, (A) drops at a lower T ~ 0.5 MeV as compared to the MBCS gap. The 
reason is that the MBCS gap incorporates the microscopic interplay between the quantal gap Aq 
and thermal one, S A. At low T the former dominates. The gap (A) takes into account only thermal 
fluctuations around the most probable value following the distribution l|94|) . The latter assumes 
equally strong coupling between A and all the intrinsic degrees of freedom, disregarding quantal 
effects. 



B. Temperature dependence of heat capacity, lelvel-density parameter, and level density 

The heat capacity C and inverse level-density parameter K — A/a obtained within the BCS and 
MBCS theories are shown in Fig. 8 as functions of T. The heat capacity usually serves as an 
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indicator for phase transitions. Within the BCS theory, a sharp discontinuity in C is seen at T = T c , 
where the gap collapses. Together with the collapse at T c of the BCS gap as an order parameter, this 
behavior of the heat capacity is a clear signature of the second-order phase transition |ll| . However, 
within the MBCS theory, this phase transition is washed out so that the temperature dependence of 
the heat capacity is a smooth curve with only a slight effect of the bending of the pairing gap in the 
region 0.5 < T < 2 MeV. A similar feature of the heat capacity has been recently reported for iron 
isotopes within the shell- model Monte-Carlo approach [17]]. Since both the order parameter A and 
the heat capacity are now continuous functions we can say that no phase transition actually occurs. 
At high T both the MBCS and BCS results approach each other. 

The inverse level-density parameter K = A/a obtained within the MBCS theory is larger than that 
obtained within the BCS theory at T < 3 MeV. At higher temperatures both theories predict almost 
the same K . Except for the low-temperature region (below T c ), where the Fermi-gas formula (|85|l is 
not valid, K increases with increasing T at T > 1 MeV, and enters the region of the experimentally 
extracted values between 8 ~ 12 MeV at T > 2.5 MeV 0. At T m 1 MeV, the value of K 
predicted by the MBCS theory is around 6 MeV, which is about twice larger than that given by the 
BCS theory. 

Shown in Fig. 9 is the logarithm of level density p(N, Z) (|88|) as a function of T. The BCS result 
shows a kink at T = T c , while the MBCS result is a smooth curve, which increases monotonously as 
T increases, exposing no signal of phase transition. At T > 2 MeV both the BCS and MBCS results 
practically coincide. 

V. CONLUSIONS 

This work has derived the modified HFB (MHFB) theory at finite temperature, which conserves 
the unitarity relation of the generalized particle-density matrix. This has been done by including 
the thermal fluctuation of quasiparticle number microscopically in the quasiparticle-density matrix. 
It has been shown that the latter can also be obtained by applying the secondary Bogoliubov 
transformation discussed in Refs. [l3L IbH]. The MHFB equations at finite temperature have been 
then derived following the standard variational procedure used in Ref. Its BCS limit yields 

the modified BCS (MBCS) equations, which have been derived previously in Refs. p^L Il4| using 
the above-mentioned secondary Bololiubov transformation. Apart from being able to restore the 
unitarity transformation, this secondary transformation helps the modified QRPA to completely 
restore the Ikeda sum rule for Fermi and Gamow- Teller transitions, which has been violated within 
the renormalized QRPA. 

The illustration of the MHFB theory has been presented within the MBCS theory by calculating 
the neutron pairing gap and thermodynamic quantities for 120 Sn. Detailed analyses of the results 
obtained show that the calculations for open-shell nuclei within the MBCS theory need to be carried 
out using the entire single-particle spectrum, which includes both bound and quasibound levels 
in a large configuration space of about 7 major shells up to 126 - 184 one. When the use of a 
reduced spectrum is unavoidable, the reduction should be done symmetrically from both side of the 
chemical potential \x so that the distribution of quasiparticle-occupation number can be properly 
taken into account as it is symmetric with respect to \x. The MBCS gap decreases monotonously with 
increasing T and does not vanish even at T ~ 5 MeV. The discontinuity in the BCS heat capacity at 
the critical temperature T c is also competely washed out, showing no signature of superfluid-normal 
phase transition. The temperature dependences of level density and level-density parameter are also 
smooth. 

The behavior of the MBCS gap as a function of T is found in qualitative agreement with that 
given by the macroscopic treatment using the Landau theory of phase transitions in the sense 
that both gaps do not collapse at the critical temperature of the BCS superfluid-normal phase 
transition. However, quantitative discrepancies between microscopic and mascroscopic approaches 
are evident. In the low-temperature region, due to the microscopic interplay between quantal and 
thermal components, the MBCS gap starts to decrease at a higher T with increasing T as compared 
to the macroscopically averaged gap (A) . At high temperatures T > 2 MeV, the MBCS gap continues 
to decrease, while (A) remains nearly constant and even start to increases with increasing T. 

The MBCS equations also shows that thermal fluctuations can induce a pairing gap even for 
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closed-shell nuclei. Results obtained using a single-particle space restricted to hole orbitals have 
shown that such a thermally induced gap increases with increasing temperature. However, its 
magnitude is negligible compared with the gap in open-shell nuclei. Therefore, it can be safely put 
to be equal to zero at T < 5 - 6 MeV. 
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APPENDIX A: RESTORATION OF IKEDA SUM RULE WITHIN THE MODIFIED 

QRPA 



At T = 0, in general, if the quasiparticle correlations are significant so that the correlated ground 
state |0) deviates appreciably from the quasiparticle vacuum or QRPA vacuum, the secondary 
Bogoliubov transformation (|47|) can be used to derive a symmetry-conserving theory, which treats 
the ground-state correlations within a microscopic and self-consistent framework. In this case the 
quasiparticle occupation number rij , which characterizes the magnitude of the ground-state correla- 
tions, can be evaluated from the renormalized QRPA backward-going amplitudes as has been 
discussed thoroughly in Refs. 0, Izl. In this section we will prove that the modified QRPA 
theory, which has been derived in Ref. |l3| using the secondary Bogoliubov transformation in the 
form of Eq. J73J), indeed conserves the Ikeda sum rule. 

1. Ikeda sum rule 



The Ikeda sum rule for Fermi (J = 0) and Gamow- Teller (J = 1) transitions is defined with 
respect to the ground state \g.s.) of the final nucleus (N, Z) as 

S~-S+ = J2 \{Ji\p-\g.*.)\ 2 - £ \(Ji\[3 + \g.s.)\ 2 = (2J+l)(N - Z) , (Al) 

i i 

where the squared /^"-transition matrix element \{Ji\f3~\g.s.)\ 2 is calculated as 



\(Ji\p-\g.8.)\ 2 = 



(Ji\ J2 ^3>3* v 3„ A hM M ) + u ^ v i^ ju {JM)]\g.s.) 



(A2) 



and P + = (/? )t. The notation Ojj { = (— ) J m Oj-m is used hereafter, and q^j v denotes single- 
particle matrix elements corresponding to the Fermi or Gamow- Teller transition. The subscripts ir 
and v denote proton and neutron, respectively. The quasiparticle-pair operators A'*- (JM) and 
A j7ijv (JM) are 



A* jwjv (JM) = (j*rn„j v m v \JM)a^ 



A jr {JM) = [A jr {JM)^ 



(A3) 



Their exact commutation relation is 



[A jUL (J'M'),Al jv (JM)] = Sjj'Smm'Sm'Jmi 



~ S 3*3'tt zJ U^ m njlm^\J'M'){j^m„j v mu\JM)al vmv aj^ m i (A4) 



~ S 3„3l (?Tr m Utymv\J'M')(j n mnj v m v \JM)a^ mn aj^ m ^ . 
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2. Fulfillment of Ikeda sum rule within QRPA 



The QRPA treats the excited state \Ji) as a one-phonon state 

\Ji)=Q i jMi\KPV > (A5) 
while the ground state |RPA) of an even-even nucleus is treated as the phonon vacuum 

QjMi|RPA) = . (A6) 
The 7rz/-phonon operator Q\ Mi is defined as 

Q\kh = E l X J J 2 A L, ( JM ) - Y Lf A M* ( > QjMi = [Q] }A J ■ (A7) 



In order to obtain a set of linear equations with respect to the amplitudes and Yj' 1 ^, the 

QRPA assumes the quasiboson approximation, which neglects the contribution of the last two terms 
~ s a at the r.h.s of Eq. I|A4(1 in the average over the ground state |RPA), i.e. 

(RPA|[A i;A (J'M'),At^(JM)]\KPA) a Sjj>S MM '6 Mir 5 jvjL . (A8) 

Within this approximation (| A8|> . the condition for the 7w-phonon operators (|A7|) to be bosons, i.e. 
satisfy the commutation relation 

(RPA\[Q JMl ,Ql, M , l ,}\RPA) = Sjj, Smm'Sw , (A9) 
leads to the following normalization relation for the amplitudes X^ J% } and Yj J jl 



Using the inverse transformation of (| AT|> . one expresses Aj^(J M) and Aj^j v (JM) in terms of the 

7w-phonon operators Q* JMi and Qjmi- Substituting the result into Eq. l|A2fl and using it to evaluate 
the left-hand side (l.h.s.) of Eq. (|A1|) we obtain 

(S- - 5 + ) RPA - E l(RPA|QjAfi/?~|RPA)| 2 - E I(RPA|Qjm 2 /3 + |RPA>| 2 



= I E Qi«jA u i« v 3 V X M« + v h u i» Y hjJ\ 2 - I E QM*( v i-« u iv X Mv + u 3* v 3» y mJ\* 



(AH) 



= E l*wJ 2 « - vl) = 2(2J+ 1)[E%,< -E°^ U U = (2^+l)(^ - Z) • 

In the above derivation the normalization condition (|A10|) and the usual BCS equation for the 
particle number are used together with the property J2j \9jdk\ 2 = 2(2 J + 1)0^ (i = (tt,^), k — 
(v, 7r)) for the single-particle matrix elements of Fermi and Gamow- Teller transitions. This derivation 
shows that the QRPA fulfills the Ikeda sum rule. 



3. Violation of Ikeda sum rule within the renormalized QRPA 



By neglecting the contribution of the two last terms at the r.h.s of Eq. i|A4(l in the ground state 
(|A6I) . the quasibson approximation (| A8|) ignores the Pauli principle between the quasiparticle-pair 
operators (|A3|) . This causes the collapse of the QRPA at a certain critical value of the interaction 
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parameter, where the solution of the QRPA equations becomes imaginary. The renormalized QRPA 
has been proposed as a method to cure this inconsistency 0, 152, |23j . 

This approach assumes that, instead of the quasiboson approximation (|A8|) . the following com- 
mutation relation holds in the average over the correlated ground state |RPA) 

(B^\[A jUl (J'M'),A]^(JM)]\^A) = s JJ' s *M' s M'jMl, D M* > ( A12 ) 

where 

= 1 - n h - n ju , «i = (RPA|«U«im|RPA) ^ . (A13) 

This means that the renormalized QRPA takes into account the contribution of the diagonal elements 
of the last two terms at the r.h.s of Eq. (|A4(I in the correlated ground state |RPA). 
The renormalized 7w-phonon operators Q JMi and Qjnn are introduced as 

QjMi = E -7^=l X %! A ljM M ) - y { /2 A J^(JM)] , Q JMi = [Q] JK J . (A14) 

The correlated ground state |RPA) is defined as the vacuum with respect to the renormalized phonon 
operators, i.e. 

QjAftlRPA) = . (A15) 

Because of Eq. QA12L these renormalized nv phonon operators satisfy the boson commutation 
relation 

(RPA|[Q JMi ,Q^ M , i ,]|RPA) = Sjj'Smm'Sw (A16) 

in the correlated ground state |RPA) provided their X and y amplitudes satisfy the same normal- 
ization condition as in the QRPA, i.e. 



E^M^-y^M^-^ . (ait) 



The phonon energies, Xj J f and y^j amplitudes are found by solving the nonlinear QRPA-likc 
equations, whose submatrices contains the factor Dj^j^ . The latter is found from the equation 

D,r i » ' V2) E^i" ^TT- + D r'f ^-} . (AM) 



ji j" 



"3 



The presence of the factor Dj w j u makes the solution of the renormalized QRPA always real as the 
interaction strength is reduced by this factor so that the collapse is avoided. However, the inverse 
transformation of l|A14|) now becomes 



4*( JM ) = VD^Yt^SiQjMi +y { j J jQjm] > Aj*jMM) = [^ > (JM)]t . (A19) 

JMi 



Using Eq. (|A19|) to evaluate the quantity S — S + in the same way as in the derivation (|Allj> . one 
finds 

(S- - s+) RRPA = \(RPA\QjMip- |RPA)I 2 - E KRPAIQjm^^+IRPA)! 2 



= E D ^,I<?^,I 2 (< 



(A20) 
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This quantity is smaller than (2 J + 1)(N — Z) since Dj 7t j u < 1. Hence, Eq. (|A20ll shows that the 
renormalized QRPA violates the Ikeda sum rule. 

We notice that, although the renormalized QRPA takes into account Eq. (|A12(1 . it neglects the 
following commutation relation between the scattering-quasiparticle pairs 

<Rr^|[flt ; ^(J'MO,S^(JM)]|RPA) = Sjj'SMM'S^ArMi* ~ n i«) > ( A21 ) 

where 

b mAJM) = - ^ m ^um v \JM)a]^a ju rni , (A22) 

The omission of the contribution of scattering-quasiparticle operators £?j, ., (J'M') and Bj % j v (JM) 
is the source that leads to the underestimation of the quantity S~ — S + within the renormalized 
QRPA. 

4. Restoration of Ikeda sum rule within the modified QRPA 

The modified QRPA makes a further step by taking into account the effects of ground-state 
correlations on the quasiparticle and collective excitations. This has been realized in Ref. |l3j | using 
the secondary Bogoliubov transformation (|73|) . where rij is the quasiparticle-occupation number in 
the new correlated ground state |RPA): 

nj = (RPA\a\ m a jm \WK) ^ . (A23) 
The modified QRPA phonon operators are introduced as 

Q\ui = E J L, ( JM ) - YjS A ^ ( J M)] > QjMi = [Q f JA J , (A24) 

where A^^(JM) and Aj^j v (JM) are the creation and destruction operators of a modified- 
quasiparticle pair 

A]^(JM) = J2 {j*m«jum v \JM)a]^a] vmv , A jf (JM) = [A jf (JM)^ . (A25) 

The new ground state |RPA) is defined as the vacuum for the modified phonon operator, i.e. 

Q|RPA) = . (A26) 

The transformation from the single-particle operators aj m and a,j m to the modified quasiparticle 
opertors aj m and Oj m is obtained after successively applying the usual (|72|l and secondary l|73|> 
Bogoliubov transformations, and has the form similar to the inverse transformation of (|72(l 

4n = *i*] m + i-y- m ^3-m , (~y- m a^ m = (-y- m uja^ m - vfa) m , (A27) 

where 

Uj = Uj yjl - n j + vj ^/Wj , Vj = Vj y/1 - n,j - u 3 . ( A28) 

Using the secondary Bogoliubov transformation (|73() to express Aj ■ (JM) and Aj v j„ (JM) in terms 
of A]^. (JM), A jlvji ,(JM), B]^(JM), and B j%ju (JM) (See Eqs. (9) and (10) of Ref. 0), we find 
that Eqs. l|A"T2"|) and \KZ\$ hold in the correlated ground state (|A"26|) if fit (JM), and B j7rJv (JM) 
commute, while ( JM), Aj^j^JM) obey the same commutation relation I|A8J| with respect 

to the correlated ground state |RPA). Therefore the set of equations to define the energy and 
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amplitudes X^^ and Y}jl of the modified 7r^-phonon excitation have the same form as of the 
usual QRPA ones. This set of equations is called the modified QRPA equations. The amplitudes 
Xf^l and obey the same normalization condition as in Eq. I|A10|) . namely 

£ {X^IMV - Yj^Y^p] = Sjj.Sh, , (A29) 

jujv 

The set of modified QRPA equations should be solve simultaneously with the normalization condition 
(|A29() and the equation for rij . The latter is evaluated from Eq. I jAlSj l as 

(A30) 

where is expressed in terms of Xj^P and Y^ as [See Eq. (16) of Ref. |l3[l] 

y$> = X^y/il - n 3 )(l n r ) + Y^y/nFi? . (A31) 

Because of the renormalization factors y(l — nj){l — rij>) and ^ynJnJT the modified QRPA can 
avoid the collapse in a way similar to that of the renormalized QRPA. Using the same secondary 
Bogoliubov transformation l)73|l , it is easy to see that the modified QRPA phonon operators contain 
all the operators A-^(JAf) and A jwju (JM), B]^(JM), and B jwjl/ (JM) in its definition [See Eq. 
(13) of Ref. [13]- This feature allows the modified QRPA to restore the Ikeda sum rule as shown 
below. 

The quantity S~ — S + is calculated within the modified QRPA as 

(S- - 5+) M rpa = ]T |(RPA|gj M4 /r|RPA)| 2 - £ |(RPA|Q JMl /3 + |RPA)| 2 

i i 



= E l^iJ 2 (4 - 4) = 2 ( 2 ^ + 1)E%4 - E^« 2 J = (2-7+ 1)(JV - Z) , 

making use of the normalization condition (|A29|1 . Eq. (|A28|I for Vj in terms of Uj and Vj, and MBCS 
Eq. (|79|l for particle number. We've just proved that the modified QRPA restores the Ikeda sum 
rule. 
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FIGURE CAPTIONS 

Fig. 1 Occupation probabilities within MBCS theory as functions of single-particle energies ej 
for neutrons at T = 0.4, 1, 3, and 5 MeV. (A thicker line corresponds to a higher T as indicated 
in panels (a) and (d)). Panel (a) shows the Bogoliubov coefficients Uj (dotted lines) and Vj (solid 
lines). Pannels (b) and (c) show the product UjVj and the difference uj — v'j, respectively. Panel 
(d) shows the quasiparticle occupation number rij. Panels (e) and (f) show the factor (1 — 2rij) and 
5Nj = y/n,j(l — n,j), respectively. The open circles marked on the lines at T = 5 MeV in (a) and 
(d) correspond to the positions of single-particle levels. 

Fig. 2 Partial quantal (Aq)j- and thermal SAj gaps as functions of single-particle energies ej 
for neutrons at T = 0.4, 1, 3, and 5 MeV. A thicker line corresponds to a higher temperature as 
indicated in panel (a). 

Fig. 3 Neutron pairing gap as a function of temperature T. The thick solid line represents the 
MBCS gap A. Its two components, the quantal gap Aq and thermal gap <5A, are shown by dashed 
and dash-dotted lines, respectively. The BCS gap is shown by the dotted line. 

Fig. 4 Neutron pairing gap as a function of temperature T. The thick solid line represents the 
MBCS gap A obtained using the entire single-particle spectrum, as in Fig. 3. The dashed line shows 
the result obtained using an N = 28 core. The thin solid line represents the result obtained using 
the same core and removing two highest levels lfc 17 / 2 an d ljii/2- The dash-dotted line is the result 
obtained using the same core and removing three highest levels, lfci7/2, ljii/2 3 and lji5/2- The 
dotted line shows the BCS gap. 

Fig. 5 Thermally induced pairing gap for protons in 120 Sn as a function of T. 

Fig. 6 Probability distribution P(A) for 120 Sn as a function of gap parameter A at different 
temperatures shown by the numbers (in MeV) next to the thick solid curves. The dashed curves in 
(b) represent the Gaussian distribution at various temperatures shown by the italic numbers 
(in MeV) next to these curves. 

Fig. 7 Neutron pairing gap for 120 Sn as a function of T within the MBCS theory and the 
macroscopic treatment following the Landau theory of phase transitions. The thick and thin lines 
show the gap A and the averaged gap (A), respectively. 

Fig. 8 Heat capacity (a) and the inverse level-density parameter K = A/a (b) as functions of 
temperature for 120 Sn. The thick solid lines denote results obtained within the MBCS theory. The 
dotted lines show the BCS results. 

Fig. 9 Logarithm of level density as a function of T for 120 Sn. The notation is as in Fig. 8. 
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